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1. Introduction 
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• Abstract. 

, Coulomb and charge-symmetry breaking effects in the ^He ground state within 

I the dressed dibaryon model developed recently for 2N and 3A^ forces are examined 

' in detail. Particular attention has been paid to the Coulomb displacement energy 

^T) ', AEc = -Eb(^H) — Esi^iic) and rms charge radii of and "^He. A new scalar 3A^ 

' force between the third nucleon and dibaryon is found to be very important for a 

correct description of the Coulomb energy and rms charge radius in ^He. In view of 
the new results for AEc obtained here, the role of the effects of charge symmetry 
Q ■ breaking in the nuclear force is discussed. 

> 
X 

The problem of accurate description of Coulomb effects in ^He in the current 3N 
approach of the Faddeev or variational type has attracted much attention for last three 
decades (see e.g. 012] and the references therein to the earlier works). The AEc problem 
dates back to the first accurate 3A^ calculations performed on the basis of the Faddeev 
equations in the mid-1960s 3 . These calculations first exhibited a hardly removable 
difference of ca. 120 keV between the theoretical prediction for AEq" — 640 keV and 
the respective experimental value AE^^ ~ 760 keV. In subsequent 35 years, numerous 
accurate 3A^ calculations were performed over the world using many approaches, but this 
puzzle was still generally unsolved. The most plausible quantitative explanation (but 
yet not free of serious questions) for the puzzle has been recently suggested by Nogga 
et al. |2I. They have observed that the difference in the singlet ^5*0 scattering lengths 
of pp (nuclear part) and nn systems (originating from the effects of charge symmetry 
breaking (CSB)) can increase the energy difference between and ^He binding energies 
and thus contribute to AEq. Using some realistic, currently accepted, app and ann 
values and other small corrections, Nogga et al. were able to virtually remove the 
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gap of 120 keV between the conventional 3A^ approaches (which neglect CSB effects) 
and the experimental value. However, this success depends crucially on the accepted a„„ 
value, which is not very reliable up to date (see the details in Section IV). For another 
admissible a„„ value, the explanation given for the gap in ref.j2] is invalid. Thus, one 
should look for another, alternative, explanation for the puzzle. 

In this paper, we give such an alternative explanation of the AEc puzzle and other 
Coulomb effects in ^He without any free parameter on the basis of the new model recently 
developed for the 2N and 3A^ forces by our joint Moscow-Tiibingen group jH E]. The 
model includes an inevitable dibaryon in the intermediate state dressed with vr, a, u, 
and p fields, together with the traditional Yukawa vr and 2tt exchanges, which describe 
the peripheral part of the A^A^ interaction. Being embedded into three- and many- 
body systems, this specific two-body mechanism generates an inevitable scalar three- 
body force induced by cr-meson exchange between the dressed dibaryon and surrounding 
nucleons jH]. 

In the preceding paper we formulated a new model for the 3A^ force of the scalar 
nature and tested it in 3A^ calculations. As was demonstrated in [H], this 3A^ force is 
so strong that can explain not only the 3A^ binding energy but also other important 
characteristics of the 3A^ system. This scalar 3A^ force is closely associated with the 
generation of an intermediate dressed dibaryon in the fundamental A^A^ interaction jHE]. 
The contribution of the above scalar 3A^ force between the dressed dibaryon and third 
nucleon to the total 3A^ binding energy is much higher than that of the conventional 3A^ 
force associated with the generation of intermediate A-isobars and two-pion exchanges 
between three nucleons (see e.g. the review paper and the references therein to 
earlier works). Thus, this scalar 3A^ force should primarily determinate the properties 
of the 3A^ nuclei, e.g., the rms radii of matter and charge distributions, the probability 
of the D state, the constants of the asymptotic normalization in the S- and D-wave 
components, etc. We also established that the complicated interplay between 2A^ and 
3A^ forces in the new model is primarily responsible for Coulomb effects in ^He and the 
Coulomb displacement energy AEc- 

This paper is devoted to the Coulomb and CSB effects in ^He, which provide an 
independent and important check for the consistency and adequacy of the new force 
model. We found that all basic Coulomb effects can be quite naturally explained in 
the framework of the new force model without any additional parameter. Thus, this 
independent check provides quite strong additional support for the force model used 
here. 

The paper is organized as follows. In Sect. II, we present some theoretical framework 
for treating Coulomb effects within the new force model using the isospin formalism. 
Section III is devoted to calculation results for the ^He ground state. The detailed 
discussion of the results obtained is presented in Sect. IV, while Sect. V incorporates 
the concluding remarks of the study. 
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2. Theoretical framework 

In this section, we present a necessary formalism for our variational calculations of 
3A^ systems within the framework of the multicomponent dressed-bag model (DBM). 
In addition, we discuss here the details concerning the inclusion of Coulomb effects 
and determination of observables in the 6qN channel. The DBM described in detail 
in jH |S] differs from the traditional OBE-type A^A^ interaction models primarily by 
the presence of non-nucleonic components in the nuclear wave function. Contrary to 
numerous hybrid quark- nucleon models (popular in the 1980s), it explicitly involves 
mesonic degrees of freedom inside the dressed bags. Hence, we consider the 3A^ system, 
where the traditional 3A^ channel is supplemented by three other channels involving the 
dressed dibaryon (dressed six-quark bag) interacting with the third nucleon. In the case 
of ^He, a large scalar force appears due to a-meson exchanges between the dibaryon 
and extra nucleon, and the additional Coulomb force arises because the bag and rest 
nucleon can have an electric charge. This new Coulomb three-body force is responsible 
for a significant part of the total ^He Coulomb energy (this three-body Coulomb force 
has been missed fully in previous 3A^ calculations within hybrid 6giV models [Hj). It 
should be emphasized here that the contribution of this three-body Coulomb force to 
the total three-body binding energy is (as will be demonstrated below) quite significant 
(~ 100 keV) and makes it possible to explain, in essence, the experimental AEc value. 

The second feature of the interaction model used here is the absence of the local 
NN short-range repulsive core. The role of this core is played by the condition of 
orthogonality to the confined 6g states forbidden in the A^A^ channel. (These states 
can be identified, e.g., with locked colour 6g states having the tetraquark-diquark 
structure jH].) This orthogonality requirement imposed on the relative-motion A^A^ 
wavefunction is responsible for the appearance of some inner nodes in this wavefunction 
(very stable under variation of the A^A^-channel energy^) and respective short-range 
loops. These short-range nodes and loops lead to numerous effects and general 
consequences for the nuclear structure (see below). One of these consequences is a rather 
strong overestimation of the Coulomb contribution when using the Coulomb interaction 
between point-hke nucleons. Thus, it is necessary to take into account the finite radius 
of the nucleon charge distribution. 

2.1. Construction of a 3N variational basis and the wavefunction of the 6qN 



Here, we give the form of the basis functions used in this work and the corresponding 
notation for the quantum numbers. The total wavefunction of the 3A^ channel, "^sn, 
can be written in the antisymmetrized basis as a sum of the three components: 



I The fact that these stationary nodes play role of the repulsive core in traditional A^-force models 
has long been established [TT)] 



component 




(1) 
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where the label (i) enumerates one (of three) possible set of the Jacobi coordinates 
(rj,pj. Every component in eq.(0) takes the form 

N-y 

^3^ = EE^'rn$?- (2) 

7 71=1 

The basis functions are constructed from Gaussian functions and the corresponding 
spin-angular and isospin factors: 

where the composite label 7*^*^ = {Xik L Sjk S tjk} represents the respective set of the 
quantum numbers for the basis functions Q: Aj is the orbital angular momentum of 
the jk pair; is the orbital angular momentum of the third nucleon (i) relative to the 
center of mass of the jk pair; L is the total orbital angular momentum of the 3A^ system; 
Sjk and tjk are the spin and isospin of the jk pair, respectively; and S is the total spin 
of the system. We omit here the total angular momentum J = 1/2 and its 2;-projection 
M, as well as the total isospin of the system T = 1/2 and its projection (in this work, 
we neglect the very small contribution of the T = 3/2 component). The spin-angular 
and isospin parts of the basis functions are taken in the form 

= \{Xik : L}{s,Sk{S,k)si : S} : JM) (4) 
Xy« = Tt^i^ = \tjtkitjk)ti : TTj) (5) 

The nonlinear parameters of the basis functions and (3^n are chosen on the 
Chebyshev grid which provides the completeness of the basis and fast convergence of 
variational calculations ^T]. As was demonstrated earlier [T2j, this few-body Gaussian 
basis is very flexible and can represent rather complicated few-body correlations. 
Therefore, it leads to quite accurate eigenvalues and eigenfunctions. The formulas 
for matrix elements of Hamiltonian (for local A^A^ interactions) with antisymmetrized 
Gaussian basis are given in paper |I3J. 

Having the three-nucleon component "^sn found in the variational calculation, one 
can construct the 6gA^-channel wavefunction ^'g^^, which depends on the coordinate 
(or momentum) of the third nucleon and the a-meson momentum and includes the bag 
wavefunction (see eq.(33) of ref. {Bj). Integrating the modulus squared of this function 
with respect to the meson momentum and inner variables of the bag, one obtains the 
density distribution of the third nucleon relative to the bag in the 6gA^ channel. This 
density can be used to calculate all observables whose operators depend on the variables 
of the nucleons and bag. However, it is much more convenient and easier to deal with the 
quasi-wavefunction of the third nucleon in the 6gA^ channel, which has been introduced 
by eq.(39) of ref. p. 

To calculate matrix elements of 3BF Coulomb and OPE forces one needs the spin- 
isospin part of 6qN components of the total wavefunction. Here we give them explicitly. 
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The potential form factors in nucleon channels ip£, ' now include the spin-isospin part 
with quantum numbers of the dressed bag: 

^t''' = <pj.LMyi^ipdii,, (6) 

where Sd and td denote the bag spin and isospin respectively, and Jj and Lj are the 
total and orbital angular momenta, respectively, referring just to the vertex form factor 
(the additional letter "i" is introduced to distinguish these quantum numbers from the 
respective total angular momentum J and orbital angular momentum L of the whole 
system). Since the present version of the DBM involves bag states with zero orbital 
angular momentum (although a more general treatment can also include dressed-bag 
components with 7^ 0), we have Jj = Sd, while the bag spin and isospin are opposite 
to each other: td + Sd = 1. The isospin part of the form factor is 

The total overlap function (with its spin-isospin part) 

xi!''i^) = {^i^'mN) (7) 



can be written, e.g., as (cf. eq.(35) in ref. [U]) 




Here, J and M are the total angular momentum of the 3A^ system and its projection. A, 
and i7 are the orbital and total angular momenta of the third (ith) nucleon respectively, 
Ti is isospinor corresponding to the third nucleon. In the present calculations of the 

ground states of and ^He (</ = 1/2), we considered only the two lowest even partial 
waves {S and D) in 3A^ wavefunctions. Therefore, A, can take only two values or 2. 
Moreover, the total angular momentum of third nucleon J' is uniquely determined by 
value of Aj: J' = 1/2 at Aj = and J' = 3/2 at Aj = 2. So, we did not include the sum 
over J' in eq.®- 

Now we redefine the quasi- wavef unction of the third nucleon in 6qN channel, 
including in it spin-isospin part of the bag: 

where | JjMj) and Tt^t^^ are spin and isospin parts of the bag function respectively. 

It is easy to see that the three form factors y?^^ used in this work (^Jq, i^l, and ipl) 
determine five radial components of the overlap function ^yj{qi) and five respective 
components of the quasi- wavefunction for the 6qN channel. To specify these components 
it is sufficient to give three quantum numbers, e.g. Sd, Aj and Lj, and we will use notation 
^Sd^i Liili) these radial component: 
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At last, we give the formula for the total quasi-wavefunction of 6qN{i) component, 
separating explicitly its spin-angular and isospin parts: 

<iv = E ■■ JM) U : TT,). (10) 

^iSd \ Li ) 

The explicit dependence of this function on the isospin projection is important for 
calculation of Coulomb matrix elements and rms charge radius. 

The interaction matrix elements include the overlap integrals of the potential form 
factors with the basis functions = where all five above components 

of the overlap function enter into the matrix elements independently (certainly, some of 
the matrix elements can vanish). The explicit formulas for the above overlap functions 
and detailed formulas for the matrix elements of all DBM interactions will be published 
elsewhere. However, when calculating both the normalization of the &qN component 
and observables, the QqN components distinguishing only by their radial parts, i.e. 
by only Lj , can be summed. Thus, only three components of the QtqN wavefunction 
(orthogonal due to their spin-angular parts) remain: 5"- wave singlet one {S^ — 0): 

^MN _ j 6gN 
^00 — ^00,0 

and two triplet {Sd = 1) ones: 

j 6qN _ j 6qN j 6qN 
^10 ~ ^10,0 "T ^10,2) 

= (11) 

The total weight of each of three 6gA^(i),(i=l,2,3) components is equal to 

n?iv = ll^off + ll^iff + ll^?f f ; ^ = 1,2,3. (12) 
Now, let us introduce the relative weights of individual 6qN components: 

||,T,6gAf||2 IkT/^^A^IlS ll\T/^9^l|2 

pSO _ II ^00 II pSl _ II^IO II pD _ 11^12 II /-,o\ 
^6qN - (j) , ^6qN " (j) , ^6qN " (j) , {^'^) 

^6qN ^QqN ^6qN 

After renormalization of the total four-component wavefunction, the total weight of all 
6qN components is equal to 

^6qN — 7f- U4J 

Here we assume that the 3A^ component of the total wavefunction, "^^^ , obtained from 
the variational calculation is normalized to unity while the total weight of the three- 
nucleon component ^3jv is equal to 

1 



3N — 



PeqN. (15) 



The total weight of the L>-wave component in ^He (and wavefunctions with 
allowance for non-nucleonic components is also changed: 

Pd = Pmi^ ~ P^qN) + PeqNPsqN (16) 
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Numerical values of all above probabilities for 6qN and 3A^ components are given 
below in Table 2. The total weight of all 6qN components P^qN in the 3N system, as was 
demonstrated in [Hj, is rather large and approaches or even exceeds 10%. Furthermore, 
taking into account the short-range character of these components, the more hard 
nucleon momentum distribution (closely associated with the first property) for these 
components, and very strong scalar three-body interaction in the 6qN channel, one can 
conclude that these non-nucleonic components are very important for the properties of 
nuclear systems . 

2.2. "Smeared" Coulomb interaction 

The Gaussian charge distribution that has the rms charge radius Tc and is normalized 
to the total charge z: Arc J pr'^dr = z can be written as 

p(r)=.(^)'^^-^ o.'' = \rl (17) 

The Coulomb potential for the interaction between such a charge distribution p(r) and 
a point- like charged particle has the well-known form 

y{R: a) = l = ^erf (i?v^) 

One can derive a similar formula for the Coulomb interaction between two charges Zi 
and Z2 with Gaussian distributions with different widths ai and 02 and rms radii r] 
and r^, respectively: 

V{R; ai, 02) = ^erf (i?v^), a = or a'^ = ^(r^ + ). (18) 

K «! + ^2 6 

In our calculations, we used the following charge radii for the nucleon and dibaryon: 
(rjp = 0.87 fm, 

{rc)6q = 0.6 fm.§ 

These values lead to the "smeared" Coulomb interactions of the form 
VSnK^) = ^erf(rv/^), = 1.005 fm; 

VStip) = irf(py^), a-^f = 0.863 fm; (19) 

§ This value is simply the rms charge radius of the six-quark bag with the parameters given in ref. 0] . 
The neutral a field of the bag changes this value only slightly. The evident diff'erence between the 
charge radii of the nucleon and dibaryon can be well understood as follows: the charge radius of the 3q 
core of the nucleon is taken usually as r^'' ~ 0.5 0.55 fm, while remaining 0.3 fm is assumed to come 
from the charge distribution of the tt"*" cloud surrounding the 3q core in the proton. In contrast, the 
meson cloud of the dibaryon in our approach is due to the neutral scalar-isoscalar a meson, so that the 
dibaryon charge distribution is characterized only by the charge radius of the bare 6q core. However in 
more complete model, one should incorporate also 7r-meson (and also p- and w-meson) dressing of the 
dibaryon. So that, the charge radius of the (pn) and {pp) dibaryons will be a bit larger than the value 
0.6 fm accepted here 
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2.3. Matrix elements of the three-body Coulomb force 

The Coulomb interaction between the charged bag and third nucleon in the 3A^ channel 
is determined by the three-particle operator with the separable kernel (see eq.(40) in 
ref.ini): 



3BF Y{i)( -n' . n. n' 



where pi is the relative momentum in the jk pair, qj is the third nucleon momentum, 
E is the total energy of the 3A^ system. The kernel function "^""'PF/',, for the point-like 

^ i 

Coulomb interaction can be taken for the OSE interaction from eq.(53) of ref.[l] with 
the substitution =^ and —g'^NN =^ e^. Variational calculations require only the 
matrix elements (m.e.) of the interaction operator between the basis functions chosen. 
It is evident that the m.e. of the operator ()2()|1 can be expressed in terms of the integrals 
of the overlap functions XL^*^'(q«) ®- 

(Coul 3BF>, ^ E Ai;,. ,0)(1..)^ / ^ jS!^'"'"^' 

where td is the isospin of the bag (we remind that 5*^ = Jj, Sd + td = 1). For brevity, 
we omitted here the labels of the basis functions. After partial-wave decomposition (cf. 
eq.©) 



E- Eo- q^/2m 
integral (PT|) reduces to a sum of integrals of the form 



E^A'"H?)n(q), (22) 



(^c'oui)' = / MlW'iq, q'WAq') g'dg/dg', (23) 



where 



V^,{q.q') = - [ Jxiqr) -jxiq'r)rMr = --^Qa(^); z = (24) 



Tx J r vr 2qq' 2qq' 

Here, Q\{z) is the Legendre function of the second kind. Now, we will replace 
the Coulomb potential 1/r between the point-like charges in eqs. ()23ll!Mjl with the 
corresponding potential between the "smeared" charges: 



(V^Coui(«))' = jMq) I ljx{qr f'\'^^ Jx{q'r)rMr 



i^'x{q')q'dqq''dq'. (25) 



It is necessary to comment the calculations of such integrals. In the momentum 
representation, the integrals include the Coulomb singularity. Thus, they must be 
carefully integrated numerically. In practice, it is much more convenient to treat them 
in the coordinate space especially in the case of the "smeared" Coulomb interaction. 
However, the presence of the propagators {E — Eq — q^ /2m)~^ in our case requires the 
use of the momentum representation from the beginning. Hence, we calculated the 
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above Coulomb integrals as follows. Taking into account that the overlap functions 
x(q) in the Gaussian basis reduce to a sum of Gaussians, we approximated the above 
propagators by a sum of few Gaussians. Then, ip\{q) entering into eq.(j221) takes the 
form 

X+n 

^^(9) = EE^AZ?'e-'^<V, (26) 

i l=\ 

where the additional degrees of q arise due to the use of an antisymmetrized basis. Now, 
the integral (^3)1 reduces to a sum of terms involving only one-dimensional integrals: 

^ J r 2 a^+i ^ (^)fc+i 2a'-'=(^ + 1)a'-'=+V2 ^ ^ 

where are the binomial coefficients. Thus, this technique reduces the whole 
calculation of the three-body Coulomb interaction matrix to completely analytical 
formulas, which considerably simplify the variational calculation. 

2.4- Rms matter and charge radii 

In the 3A^ channel, the rms radii of the proton (r^) and neutron {vn) distributions are 
defined by the standard way: 

3 1 + r^^^ 

r^ = |^(^3iv|^^r,2|^3iv), (28) 

where ri = (2/3)pi is the distance between particle (1) and the system center of mass, 
Np = 3/2 + Tz and Nn = 3/2 — are the numbers of protons and neutrons, respectively. 
Np is equal to the total charge of the system Z. Then, the rms matter radius is equal 
to 

rl = '\ " = (^3^|r?|*3^) (29) 
The rms charge radius in the 3A^ sector is also defined conventionally: 



{rluhN = rl + Rl-^K. (30) 



where R^ = 0.7569 fm^ and R^ = —0.1161 fm^ are the squared charge radii of the 
proton and neutron, respectively. 

Further, we define the rms charge radius in 6qN channels as 



1 / , . 1 ^ t(^) 1 - t(^) 

{r!,h,N = (*^;^v -^{rl + Rl) + '-^Rl 



7 



2 



+ (1 + ii)rl + E ^t,tMt, t.) , (31) 
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where 



ri = ■ = api and = ■ = -(1 - a)p^ 

md + rriN rud + mN 



are the coordinates of the third nucleon (with number 1) and the bag in the c.m.s., is 
the operator of the third component of the bag isospin, Vt^t^ is the projector into isospin 
state of the bag with definite values of its isospin t and z-projection t^. The Rdit.tz) 
is the value of charge radius of the bag in the specific isospin state. These values, in 
general, are different, but their difference should be rather small and, thus, can be safely 
ignored in subsequent calculations keeping in mind the relatively low probability of the 
bag-component. Thus we take value Rd = 0.6 fm for the mean charge radius of the bag 
in all isospin states, except the state with tz = —1 corresponding to (nn)-bag. For latter 
(uncharged) state we put i?^(l, — 1) = 0. 

The total rms charge radius of the 3A^ system with allowance for both the three- 
nucleon and dibaryon-nucleon components thus takes the form 

(rl) = (1 - P,,N){rlhN + P6,N{rl),,N. (32) 

Similarly to eq. ()28|1 . one can define the rms radius of the proton (neutron) 
distribution in the 6qN channel as 

where the denominator determines the average number of protons N^'^^ (or neutrons 
N^i^) in the 6qN channel: 

^\:} - p{i) • ^"^^^ 

One can note that, if to neglect the difference between isosinglet and isotriplet 
wavefunctions in the 6qN channel, then these numbers are equal Np = 1/3 and 2/3 
(A^„ = 2/3 and 1/3) for and ^He respectively. 

The total rms radii of the nucleon distributions including both the 3N and three 
6qN components are 

= (1 — -P6gAf)^?P j + -P6gAf(^?P|)6gAf- (35) 

2. 5. Role of the pn mass difference 

To accomplish calculation for the accurate Coulomb displacement energy AEq = 
_E'b(^H) — Esi^iie), one should take into consideration some tiny effects associated 
with the mass difference between the proton and neutron. It is well known [2j that the 
above mass difference makes rather small contribution to the difference between ^He 
and binding energies. Therefore, it is usually taken into account in the perturbation 
approach. However, since the average kinetic energy in our case is twice the energy in 
conventional force models, this correction is expected to be also much larger in our case. 
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Hence, we present here the exphcit exact evaluation for such a correction term without 
usage of the perturbation theory. 

It should be added that in our model involving various 6qN components the similar 
effect originated from the mass difference of dibaryons in different charge states with 
Z = 0,1 and 2 should also be taken into account. This mass difference is equal about 
AM^ ~ 3 MeV. || The latter effect seems to yield a negligible correction to the AEq 
value, because the total probability of all 6qN components does not exceed 10 - 11 %, 
while the nucleon mass difference is half the AM^ value, i.e. ca. 1.5 MeV, at the 
probability of the 3A^ channel ca. 90%. 

In the conventional isospin formalism, one can consider that the and ^He nuclei 

consist of the equal-mass nucleons: 

mp + rUn 

m = — , 

2 

so that rrip = m + Am/2, rUn = m — Am/2, where Am = mp — mn- The simplest way 

to include the correction due to the mass difference Am is to assume that all particles 

in have the average mass 

2m„ + mp 1 

mn = = m Am, 

3 6 

while they have the different average mass 

2mp + mn 1 . 

rriHe = = m + - Am 

6 o 

in '^He. In spite of smallness of parameter Am/m, the perturbation theory in this 
parameter does not work. So we used the average mass rhn in calculation of and 
rfiHe in calculation of ^He. The results of these corrections are given in fifth row of 
Table 2. 



3. Results of calculations 



Here, we present the results of the 3N bound-state calculations based on two variants 
of the DBM. 

(i) In first version developed in the dressed-bag propagator includes three loops, 
two of them are of the type shown in Fig. 2 of ref. jl], in which each loop was found 
with the ^Pq model. The third loop consists of two such vertices and a convolution 
of the (T-meson and 6g-bag propagators (see the Fig. 2 in ref j^). 

(ii) In the second version, we replaced two above loops with the effective Gaussian form 
factor B{k), which describes the direct NN — 6g + a transition, i.e., the direct 
transition from the A^A^ channel to the dressed-dibaryon channel. 

Both versions have been fitted to the A^A^ phase shifts in low partial waves up to 
an energy of 1 GeV with almost the same quality. Therefore, they can be considered 
on equal footing. However, version (ii) has one important advantage. Here, the energy 

II The mass difference between baryons with different ST values is already included in our force model. 
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dependence arising from the convolution of the two propagators involved into the loop, 
i.e., the propagators of the a-meson and bare dibaryon, describes (with no further 
correction) just the energy dependence of the effective strength of the NN potential 
\^'^\E), which is thereby taken directly from the above loop integral. In contrast, in 
the first version of the model, two additional qqirTca loops give a rather singular three- 
dimensional integral for X^^\E), where the energy dependence at higher energies should 
be corrected by a linear term. The main difference between the results for both versions 
is that the energy dependence of X{E) for the second version is much weaker than that 
for the first variant. In addition, this energy dependence leads both to the decrease in 
the contribution of the 6qN component to all 3A^ observables and thus to the increase 
in the contribution of the two-body force as compared to the three-body force. 

Table 1 presents the calculation results for the two above versions for the following 
characteristics: - the weights of the 6qN channels and D wave in the total 3N function, 
as well as the weight of the mixed-symmetry S' component (only for the 3 channel) ; - 
the averaged individual contributions from the kinetic energy T, two-body interactions 
yC^N) plug ^Jie kinetic energy T and three-body force (V"*^^^^) due to one-sigma and 
two-sigma exchanges to the total Hamiltonian expectation. 



Table 1. Results of the 3N calculations with two- and three-body forces for two 
variants of the dressed-bag model. 





E 
MeV 


Pd% 


Ps' % 


PeqN % 


Individual contributions to H, MeV 


T 


T + V^(2iv) 


\/(3A^) 


3 


H 


DBM(I) g = 9.577*) 


-8.482 


6.87 


0.67 


10.99 


112.8 


-1.33 


-7.15 


DBM(II) g = 8.673*) 


-8.481 


7.08 


0.68 


7.39 


112.4 


-3.79 


-4.69 


AV18+UIX1) 


-8.48 


9.3 


1.05 




51.4 


-7.27 


-1.19 


3He 


DBM(I) 


-7.772 


6.85 


0.74 


10.80 


110.2 


-0.90 


-6.88 


DBM(II) 


-7.789 


7.06 


0.75 


7.26 


109.9 


-3.28 


-4.51 


AVIS+UIX^) 


-7.76 


9.25 


1.24 




50.6 


-6.54 


-1.17 



*) These values of aNN coupling constant in '^H calculations have been chosen to reproduce 
the exact binding energy of nucleus. The calculations for ^He have been carried out without 
any free parameters. 

The values are taken from jl5j . 



To compare with the respective results for the conventional A^A^ potential models. 
Table 1 also presents the results of recent calculations with the Argonne potential AV18 
and Urbanna- Illinois three-body force UIX ^3]. The Coulomb displacement energies 
AEc, together with the individual contributions to the Ai?^- value, are presented in 
Table 2. The rms radii of the charge and proton distributions in '^H and ^He found in 
the impulse approximation, as well as the respective experimental values and results 
obtained for AV18 -|- UIX A^A^ forces, are presented in Table 3. To demonstrate 
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Table 2. Contribution of various terms (in kcV) of the interaction to the '^H 
mass difference. 



Contribution 


DBM(I) 


DBM(II) 


AV18+UIX 


point Coulomb 3N only 


598 


630 


677 


point Coulomb 3N+6qN 


840 


782 




smeared Coulomb 3N only 


547 


579 


648 


smeared Coulomb 3N+6qN 


710 


692 




np mass difference 


46 


45 


14 


nuclear CSB^^ 








+65 


magnetic moments & 
spin-orbit^) 


17 


17 


17 


Total 


773 


754 


754 



1) See Table 4. 

Here we use the value for this correction from ref. [2] 



Table 3. Rms proton, Tp, and charge, rch, radii (in fm) in DBM approach 



model 


3H 


3He 






rp 


Teh 


DBM(I) 


3N 


1.625 


1.779 


1.805 


1.989 


6qN 


1.608 


1.188 


1.854 


1.412 


total 


1.625 


1.724 


1.807 


1.935 


DBM(II) 


3N 


1.613 


1.769 


1.795 


1.980 


6qN 


1.573 


1.171 


1.829 


1.396 


total 


1.613 


1.732 


1.796 


1.944 


AV18 + UIX 


1.59*) 




1.76*) 




Experiment 


1.60*) 


1.755 


1.77*) 


1.95 



*) These values are taken from ^5]. The "experimental" values of point proton radii rp have 
been obtain there from charge radii by removing the proton and neutron charge radii 0.743 fm^ 
and -0.116 fm^ respectively. 

the separate contributions of the three-nucleon and dibaryon-nucleon channels to these 
observables, we also present the values calculated separately with only nucleonic and 
6qN parts of the total wavefunction. 

We present here a few comments concerning the results in Tables 1-3. 

Comments to Table 1. 

(i) It is seen that an admixture of the mixed-symmetry 5" component in our 3iV 
wavefunction is almost half that for the conventional force model (e.g., AV18+UIX). 
This difference can be attributed to the fact that the relative contribution of two- 
body interactions to the total 3A^ binding energy in our approach is much less than 
that in the conventional force models (this follows from the results presented in the 
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seventh and eighth columns of Table 1). While it is well known that the weight of 
the <S" component is proportional to the difference between two-body spin-singlet 
and spin-triplet A^A^ interactions, the leading contribution of three-body forces in 
our approach comes from the scalar-isoscalar 3BF that is completely insensitive to 
the above difference. Therefore, as a result of this redistribution of various force 
components, the weight of the S' component decreases by almost half. A similar 
but weaker diminution with respect to conventional force models is also seen in the 
weight of the D-wave component Pd- This decrease has the same origin, viz. the 
suppression of the two-body force contribution and the large increase in the scalar 
three-body force contribution. 

(ii) Another remarkable distinction from the conventional force models is the large 
increase in the average kinetic (and potential) energy, viz. 112 keV vs. 50 keV in 
conventional force models (see the sixth column in Table 1). This increase is caused 
by the appearance of the short-range radial nodes and respective loops in the radial 
3iV wavefunctions (see Fig. 5 in the preceding paper). This large enhancement in 
nucleon velocities will strengthen all the effects associated with the nucleon currents, 
relativistic effects, meson-exchange contributions to electromagnetic observables, 
etc. 

Comments to Table 2. Here, we emphasize three important points. 

(iii) First, it is seen quite a large contribution from the Coulomb three-body force (cf. 
the differences between the entries "Coulomb 3iV only" and " Coulomb 3iV -|- QqN" 
in this table). The second and third rows correspond to the Coulomb interaction 
between point-like charges, while the fourth and fifth rows include results for the 
Coulomb interaction between properly smeared charge distributions. In both cases, 
the contribution of the three-body Coulomb interaction (which has been completely 
overlooked in previous works) is as large as ca. 110 - 240 keV and, along with other 
minor effects, can quantitatively explain the Coulomb displacement energy of ^He. 

(iv) The second point, which is closely interrelated to the first one, is rather high 
sensitivity of all above Coulomb contributions to the smearing of the charges (both 
for the proton and 6g dibaryon) with the Gaussian distribution. Table 2 shows 
that the inclusion of the smeared charge distribution reduces the Coulomb two- 
nucleon force contribution by 51 keV for both versions of the model. It should be 
compared to a difference of 29 keV in the Coulomb interactions between point-like 
and smeared nucleons for the AV18 -|- UIX force model. Smearing the charges leads 
also to significant reducing of the three-body Coulomb contribution: from 242 to 
163 keV for version I and from 152 to 113 keV for version II. However, even in 
the minimum scenario, we obtain an additional contribution of 113 keV from the 
three-body (smeared) Coulomb force. 

(v) The third interesting feature, which is distinguished from the conventional model 
result, is a quite large effect of the (small) up mass difference on the 3A^ Coulomb 
displacement energy. This effect is about twice the respective contribution for the 
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AV18 + UIX force model. This enhancement is attributed to the much increased 
average kinetic energy in our approach. Thus, the variation of this energy due to 
the np mass difference should be also much larger. 

We also add a small correction due to electromagnetic interactions and spin-orbit 
electromagnetic interaction (as in conventional models, we take a value of 17 keV 
for this correction). Including all these corrections, we obtain the total value 
AE^j = 773 keV for version I and AE^ = 754 keV for version II. Thus, we found 
a quite small space for nuclear CSB effects: -9 keV for DMB(I) and +10 keV for 
DBM(II). These values should be compared to a significant value of 65 keV for 
the AV18 + UIX force model. The more detailed discussion of CSB effects in our 
approach in the next section further corroborates this important conclusion: the 
admissible value of CSB effects in the DBM is noticeably smaller than that in the 
conventional force models. 

Comments to Table 3. 

(vi) The rms radii of the proton and charge distributions are presented in Table 3 for two 
versions of our model in comparison with the results for the A18 + UIX force model. 
As is seen in Table 3, the rms charge radii for the 6qN component in both and 
^He are much smaller than those for the 3A^ component (as could be expected in 
advance). On the other hand, the rms charge radii for the 3A^ component turn out 
to be larger than the respective experimental values in both and ^He. Thus, it is 
the contribution of the 6qN component to the total wavefunction of the 3N system 
that provides quite good agreement of the rms charge radii with the respective 
experimental values. 

A small underestimation of charge radii (especially for ^H) can be due to too small 
value for charge radius of the bag (0.6 fm) accepted in our calculations. This value 
is, in fact, the quark-core radius of the bag, but our estimate shows that the pion 
cloud will increase its charge radius up to 0.65 -r- 0.68 fm. This will lead to increase 
of charge radius of QqN component and, therefore, to some increase of the total 
charge radius. Besides that, there is a contribution to charge radius from the model- 
dependent two-nucleon current operator. As it shown in ref . , this contribution 
for AV18-FUIX force model is about 0.014 fm for and 0.009 fm for ^He. 

4. Discussion 

Here, we will discuss the main results found in the work in the general context of 
few-body physics and compare them with the respective results based on conventional 
force models. Particular attention will be paid to some general conclusions that can be 
derived from the results presented here. Let us begin with the results for the Coulomb 
displacement energy AEc = -Eb(^H) — Esi^iie). We emphasize three important points, 
where our results differ from those for conventional models. 

(i) First, we found a serious difference between conventional and our approaches in the 
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short-range behaviour of wavefunctions even for the nucleon channeL Conventional 
3A^ wavefunctions are strongly suppressed along all three interparticle coordinates 
rij due to the short-range local repulsive core, while our wavefunctions (in the 
3A^ channel) have stationary nodes and short-range loops along both all and 
the third Jacobi coordinates pk- Such a node along the p coordinate presents 
also in the 6qN relative-motion wavefunction (see Fig. 5 in ref. jH]). This very 
peculiar short-range behaviour of our wavefunctions leads to a strong enhancement 
of the high-momentum components of nuclear wavefunctions, which is indicated by 
various modern experiments, e.g., ^B.e{e,e'pp) ^Tj or pp pp'y etc. where high 
momentum transfers appear. On the other hand, these short-range radial loops 
lead to significant errors when using the Coulomb interaction between point-like 
particles within our approach. Hence, we must take into account the finite radii of 
charge distributions in the proton and 6q bag. Otherwise, all Coulomb energies are 
overestimated. 

(ii) Another important effect following from our calculations is a quite significant 
contribution of the 6qN component to AEc- In fact, just this interaction, 
which is completely missing in conventional nuclear force models, makes the main 
contribution (ca. 100 keV) to filling the gap in AEc between conventional 3A^ 
calculations and experiment. 

The large magnitude of this 3A^ Coulomb force contribution is explained by two 
factors: first, a rather short average distance (p^) between the 6q bag and third 
nucleon (which enhances the Coulomb interaction in the 6qN phase) and, second, a 
significant weight of the 6qN components where the bag has the charge -f-l (i.e., it is 
constructed from an np pair). This specific Coulomb repulsion in the 6qN channel 
should appear also in all other nuclei where the total weight of such components 
is about 10% and higher. Therefore, it should strongly contribute to the Coulomb 
displacement energies over the entire periodic table and could somehow explain the 
long-term Nollen-Schiffer paradox |TH] in this way. 

(iii) The third specific effect that has been found in this study and contributes to the 
quantitative explanation of AEc is a strong increase in the average kinetic energy 
(T) of the system. This increase in (T) has been already discovered in the first early 
3A'^ calculations with the Moscow A^A^ potential model jTH] and results in a similar 
nodal wavefunction behaviour along all interparticle coordinates but without any 
non-nucleonic component. 

The increase in (T) leads to the proportional increase in the np mass difference 
correction to AEc- As is seen in Table 2, this correction in our case is not very small 
and contributes significantly to AEc- Many other effects attributed to increasing 
the average kinetic energy of the system will arise in our approach, e.g., numerous 
effects associated with the Fermi motion of nucleons in nuclei. 
It is worth also to estimate here the possible contribution of 9g-bag component in 
^He to the AEc, which has been omitted in our present calculations. If one assumes 
that the radius of the 9g-bag is near to that for 6g-system, i.e. rgq ~ 0.6 fm, the 
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Coulomb energy of this bag (with the charge Z = 2) Ze'^/rgg is about 3 MeV. The 
probabihty of the 9g-bag is expected to be around 0.1 - 0.2 %. So that the additional 
Coulomb energy contribution to AEc should be about 3-(0.lH-0.2)-10-3 ~ 3^6 keV, 
i.e. a quite small value as compared to the quantity 120 keV coming from 6qN 
component. It can also be expected that the corrections to kinetic energy difference 
due to contribution of 9q component would be very small. 

The best explanation for the AEc value in the framework of conventional force 
models published up to date j2! is based on the introduction of some CSB effect, i.e., 
the difference between nn and pp strong interactions. At present, two alternative values 
of the nn scattering length are assumed: 

= _i8.7 fm, and 4n = -16.3 fm. (36) 

The first value has been extracted from the previous analysis of experiments 
d{7i~,'j)nn ^3 (see also ref. ^efs. therein) and is used in all current NN 

potential models, while the second value in (j^Bjl has been derived from numerous three- 
body breakup experiments n + d ^ nnp done for the last three decades. In recent years, 
such breakup experiments are usually treated in the complete Faddeev formalism, which 
includes most accurately both two-body and three-body forces [22] • Thus, this a„„ value 
is considered as quite reliable. However, the quantitative explanation for the AEc value 
in conventional force models uses just the first value of essential point of all 

the construction. At the same time, the use of the second value a„„(= —16.3 fm) (which 
is not less reliable than the first one) fails completely the above explanation! 

Therefore, in order to understand the situation more deeply and to determine the 
degree of sensitivity of our prediction for AEc to variation in a„„, we made also our 
3A^ calculations with two possible values of a„„ from eq. ()36|) . These exact calculations 
have been carried out with the effective values of the singlet-channel coupling constant 
corresponding to the V^gN part of the A^A^ force: 

ff 1 1 2 

"^SRei ^o) = 2 + -Xnp] (37) 

ff 1 1 2 

Xsjji So) = -Xnn + --^rap- (38) 

In the above calculations, we use the value A„p = 328.9 MeV that provides the accurate 
description of the ^5*0 np phase shifts and the experimental value of the np scattering 
length a„p = —23.74 fm Here, we employ the value App = 325.523 MeV fitted to the 
well-known experimental magnitude app = —8.72 fm and two A„„ values corresponding 
to two available alternative values of the nn scattering length: anl = —16.3 fm and 

(2) 

ann = —18.9 fm. The calculation results are presented in Table 4. 

As is seen in Tables 2 and 4, within the DBM (version I) one has AEc= 773-18=754 
keV, so that the version I of DBM can reasonably reproduce the the experimental 
Coulomb displacement energy AEc with the lower (in modulus) value a„„ = —16.3 fm, 
while this model overestimates AEc by 54 keV with the larger (in modulus) value 
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Table 4. Contribution of charge symmetry breaking effects to the '^H 
difference. 





AEc, keV 


DBM(I) 


DBM(II) 


-16.3 


-18 


-39 


-18.9 


+45 


+26 



o-nn = —18.9 fm. Thus, the DBM approach, in contrast to the conventional force models, 
prefers the lower (in modulus) possible value -16.3 fm of the nn scattering length. 

Now, let us discuss shortly the magnitude of CSB effects in our model. The 
difference between a„„ and so-called "pure nuclear" pp scattering length is usually 
considered as the measure of CSB effects at low energies. The value is extracted from 
pp scattering data when the Coulomb potential is disregarded. The model dependence 
of the latter quantity was actively discussed in the 1980s 1^ OH] . However, the 
majority of modern NN potentials fitted to the experimental value a^p = —8.72 fm give 
the value = —17.3 fm when the Coulomb interaction is discarded. It is the value 
that is adopted now as an "empirical" value of the pp scattering length [22] ■ Thus, 
the difference between this value and a^n is usually considered as the measure of CSB 
effects. However, our model (also fitted to the experimental value app = —8.72 fm) gives 
a quite surprising result for pp-scattering length when the Coulomb effects are removed: 

a^p{DBM) = -16.57 fm, (39) 

which differs significantly from the above conventional value (by 0.8 fm) due to the 
explicit energy dependence of the A^A^ force in our approach. 

Thus, if the difference — ann is still taken as the measure of CSB effects, the 
smallness of this difference obtained in our model testifies to a small magnitude of the 
CSB effects, which is remarkably smaller than the values derived from conventional OBE 
models for the NN force. 

Now, let us pass to the data from Table 3 for the radii of the charge and proton 
distribution in ^H and ^He. It is seen that both our versions (DBM(I) and DBM(II)) 
give quite similar values for all radii. The most interesting point here is the importance 
of 6qN component contributions. In fact, the contribution of the 6qN channel shifts all 
radii, i.e., Vch and rp in ^H and ^He, predicted with pure nucleonic components, much 
closer to the respective experimental values. 

Thus, the dibaryon-nucleon component also works in a right way in this aspect. It 
is interesting to note that, in general, the predictions of our two-phase model are quite 
close to those of the conventional single-phase AV18 + UIX model. This means that 
(at least for many static characteristics) our multi-channel model is effectively similar 
to a conventional purely nucleonic model. However, this similarity will surely hold 
only for the characteristics that are sensitive mainly to low momentum transfers, while 
the properties and processes involving high momentum transfers will be treated in two 
alternative approaches in completely different ways. 
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'He mass 
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It is worth here to add a few remarks about role of the Coulomb effects in 3A^- and 
4A^-continuum and its interpretation on the basis of the present model. There are a 
few long-standing puzzles in the field which still cannot be resolved within conventional 
models of 2N and 3A^ forces. Here are: 

(i) The values of Ay and Tn in nd and pd radiation capture at low energies are strongly 
underpredicted by conventional theoretical approach fI7\ . 

(ii) The pronounced discrepancy for the fore-aft asymmetries in mirror reactions 
3H(7,n)2H and 3He(7,p)2H at energies a few MeV above the thresholds and also 
for inverse capture processes fl^ EH] • 

(iii) The ratio of cross sections (7,p)-to-(7, n) for ^He at energies a few MeV above the 
thresholds is as large as 1.7 while the conventional theory predicts only the values 
ca. 1.3 - 1.4 ISO]. 

The discrepancy in (i) can be reduced somehow by inclusion of the strong two- 
body meson-exchange currents incorporating A-current with (unnaturally) high cut-off 
parameter A^ata and Aj^tva jSll ISSj . The puzzle in (iii) could be explained or reduced by 
assuming a strong charge-symmetry breaking force component which is in an evident 
contradiction with results of other experiments and also with conclusions of the present 
force model. 

It should be stressed that all three above-mentioned puzzles are interrelated to 
contribution of P-states in nd and especially nT systems. For example, as has been found 
in recent 4N calculations the P-wave peak in uT elastic scattering at energies ca. 3 - 
4 MeV cannot be explained by the fully realistic 4N-calculations within the conventional 
force model. In addition, the issue (iii) can be explained by an enhancement of the 
Coulomb effects in p—^B. exit channel. 

As follows from the results of the present and preceding works, our force model 
predicts inevitably an enhancement of the P-wave contributions in Nd and especially 
NT near continuum (due to additional strong scalar 3A^-force) and also an enhancement 
of the Coulomb effects in pd and p^H near continuum. So the present approach 
could remove or reduce noticeably the above discrepancies in 3A^ and 4A^ low-energy 
continuum. 

5. Conclusion 

Here, we will summarize the main results of this work. In the previous work, we fixed 
the only coupling constant, QaNN, to obtain the experimental value of the triton binding 
energy. Then, all other calculations in both previous and this works did not include 
any fitting parameter. Thus, their results can be considered as a stringent test for the 
proposed new model for 2N and 3N forces. 

First, we point to the precise value obtained for the Coulomb displacement energy 
AEc of the A = 3 system in the developed model. It should be emphasized that, 
contrary to other studies based on conventional force models (using the 2N and 3N 
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forces generated via the meson-exchange mechanism), this explanation does not require 
any noticeable CSB effect, ahhough our model is still compatible with such effects. 
However, these CSB effects do not contribute remarkably to ^Eq in our approach. 
Two basic sources of this contribution, which differ from conventional force models, 
should be indicated here: 

- the three-body Coulomb energy of the interaction between the dressed bag and 
third nucleon; and 

- quite significant correction to the kinetic energy of the system due to the np mass 
difference and high average kinetic energy. 

The second general point that must be emphasized is a rather large admixture 
of dibaryon-nucleon components in both and '^He, which has been calculated in 
a completely consistent way. Closely associated with the above QqN components, it 
is a specific energy dependence of the two-body force in a three- (and many-) body 
system. This energy dependence strongly reduces the contribution of two- body force 
when a strong attractive three-body force is added to the system Hamiltonian. This is a 
manifestation of a very specific new interplay between two- and three-body forces: the 
stronger the three-body force, the smaller the total contribution of the two-body force to 
the nuclear binding energy! By this way, a very natural density dependence of nuclear 
interactions appears from the beginning. Thus, the general properties of the 3A^ system, 
where forces so much differ from any conventional model force, should appear also much 
differ from the predictions of any conventional model and, hence, from experiment. 

It was very surprising to find that the characteristics of the SA/" system in 
our case turned out to be very close to the predictions of the modern force model 
(such as AV18 -I- UIX) and thus to experiment. This gives us a good test of the 
self-consistency and accuracy of the new force model. However, predictions of the 
present 2N- and 3A^-force model in other aspects will strongly deviate from those for 
conventional models. First, these are the properties determined by the high-momentum 
component of nuclear wavef unctions. The point is that the system described by our 
multi-component wavefunctions explicitly including dibaryon components can easily 
absorb quite high momentum transfers, which can hardly be absorbed by the system 
described by traditional multi-nucleon wavefunctions. Therefore, to fit the experimental 
data corresponding to large momentum transfers (~ 1 GeV/c), many types of meson- 
exchange and isobar currents are often introduced to theoretical frameworks. However, 
these currents are often unrelated to the underlying force model. Hence, it is rather 
difficult to check the self-consistency of such calculations, e.g., the validity of gauge 
invariance etc. 

Thus, the alternative description given here by the new force model can be more 
self-consistent and straightforward. One aspect of this new picture is evident - the 
present model applied to any electromagnetic process on nuclei automatically leads to 
a consistent whole picture of the process: single-nucleon currents at low momentum 
transfers, meson-exchange currents (including new meson currents) at intermediate 
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momentum transfers, and quark counting rules at very high momentum transfers, 
because the model wavefunction explicitly includes multinucleon, meson-exchange and 
multiquark components. 
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